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If, as recently proposed by Farrar and Piran, Cen A is the source of cosmic rays detected above the 
Greisen-Zatsepin-Kuz'min cutoff, neutrons are ~ 140 more probable than protons to be observed 
along its line of sight. This is because the proton flux is rendered nearly isotropic by 0{fiG) 
intergalactic magnetic fields. With the anticipated aperture of the Southern Auger Observatory, 
one may expect on the order of 2 neutron events/year above 10^^ eV in the line of sight of Cen A. 
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The energy spectrum of cosmic rays (CRs) is well fitted 
by power laws with increasing index for energies above 
4 X 10^^ eV (the "knee") flattening again above 5 x 10^^ 
eV (the "ankle" ), yielding the overall shape of a leg. Over 
the last third of the century, ingenious installations with 
large effective areas and long exposure times — needed to 
overcome the steep falling flux — have raised the tail of 
the spectrum up to an energy of 3 x 10^° eV, with no 
evidence that the highest energy recorded thus far is Na- 
ture's upper limit [|| . The origin of these extraordinarily 
energetic particles continues to present a major enigma 
to high energy physics ||] . 

The main problem posed by the detection of CRs of 
such energy (if nucleons, gammas, and/or nuclei) is en- 
ergy degradation through inelastic collisions with the uni- 
versal radiation fields permeating the universe. There- 
fore, if the CR sources are all at cosmological distances, 
the observed spectrum must virtually end with the 
Greisen-Zatsepin-Kuz'min (GZK) cutoff at £" «:! 8 x lO^^ 
eV Q . The spectral cutoff is less sharp for nearby sources 
(within 50 Mpc or so) . The arrival directions of the trans- 
GZK events are distributed widely over the sky, with 
no plausible counterparts (such as sources in the Galac- 
tic Plane or in the Local Supercluster). Furthermore, 
the data are consistent with an isotropic distribution of 
sources in sharp constrast to the anisotropic distribution 
of light within 50 Mpc [Q. The difficulties encountered 
by conventional acceleration mechanisms in accelerating 
particles to the highest observed energies have motivated 
suggestions that the underlying production mechanism 
could be of non-acceleration nature. Namely, charged 
and neutral primaries, mainly light mesons (pions) to- 
gether with a small fraction (3%) of nucleons, might be 
produced at extremely high energy by decay of super- 
massive elementary X particles {mx ^ 10^^ — 10^* eV) 
1^. However, if this were the case, the observed spectrum 
should be dominated by gamma rays and neutrinos, in 
contrast to current observation ^ ! Alternative explana- 
tions involve undiscovered neutral hadrons with masses 
above a few GeV neutrinos producing nucleons and 
photons via resonant Z-production with the relic neu- 
trino background Q , or else neutrinos attaining cross sec- 
tions in the millibarn range above the electroweak scale 
iPl. A controversial correlation between the arrival di- 



rection of CRs above 10^° eV and high redshift compact 
radio quasars seems to support these scenarios pO[ |. 

Over the last few years, it has become evident that 
the observed near-isotropy of arrival directions can be 
easily explained if even the highest energy cosmic rays 
propagate diffusively, camouflaging a unique source only 
a few Mpc away Within this framework, the par- 

ticles experience large deflections through randomly ori- 
ented patches of strong magnetic fields 0{fj,G) p2| , p^ . 
Recently, Farrar and Piran (FP) |l4| noted that an extra- 
galactic magnetic field of ~ 0.3 fiG would bend CR paths 
sufhciently, allowing one to trace back trans-GZK orbits 
from the Earth's northern-hemisphere to the southern 
radio galaxy Cen A. Moreover, they show that the flux 
of Cen A at 10^^ eV (at Earth) is comparable to that of 
all other sources in the universe, and assuming a diffuse 
propagation of particles above this energy they predict 
a CR anisotropy of order 7% (or less). Both estimates 
strongly support the single-source hypothesis. If this is 
the case, and the absence of the GZK cutoff is a reflection 
of our coincidental position near Cen A (d w 3.4 Mpc), 
it must be that the emission of uncharged particles from 
Cen A should render an enhancement of the CR flux in 
the southern hemisphere. 

Cen A is a complex Fanaroff-Riley (FR) I radio- 
loud source {I « 310°, b « 20°) identified at optical fre- 
quencies with the galaxy NGC 5128 ||l^. The radio mor- 
phology is intricate with large non-thermal radio lobes. 
In particular, the structure of the northern middle lobe 
resembles the "hot spots" which exist at the extremities 
of FR-II galaxies , although for Cen A the brightness 
contrast (hot spot to lobe) is not as extreme as in e.g. 
Cyg A |l8|. The energetics of acceleration in hot spots 
were discussed in ||l9| . The criteria were applied in |po|| to 
show the plausibility of attaining trans-GZK energies in 
the hot spot of Cen A. Moreover, EGRET measurements 
1^ of the gamma ray flux for energies > 100 MeV allow 
an estimate ~ 10^"'^ erg for the source [^. This 
value of Lj is consistent with an earlier observation in 
the TeV-range during a period of elevated activity , 
and is considerably smaller than the estimated bolomet- 
ric luminosity Lboi ~ lO'^^erg s~^ [ p^ . 

CR "lore" convinces us that the TeV 7-ray emission 
is a result of synchrotron radiation of electrons or pro- 
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tons of still higher energy ||2J,|2^. Strictly speaking, the 
observed 7-radiation is related to: (i) the development 
of pairs cascades triggered by secondary photopion prod- 
ucts that cool instantaneously via synchrotron radiation 
(ii) the synchrotron radiation of protons itself that be- 
comes a very effective channel to produce high energy 
7-rays above 10^^ eV. There are plausible physical argu- 
ments 11^:^ as well as some observational reasons ^] 
to believe that when proton acceleration is being limited 
by energy losses, the CR luminosity LcR ~ L^. The low 
ratio Lj/L\jo\ thus leads us to assume that both ultra 
high energy CR and 7 production take place in the lobes 
with the bulk of the softer radiation coming from the 
core. 

Following FP we introduce e, the efficiency of ultra 
high energy CR production compared to high energy 7 
production — from the above, we expect e ~ 1. Using 
equal power per decade over the interval 1 x lO^^eV < 
E < 4 X lO^^eV, we estimate a source luminosity 



dEdt 



1.7 6^41 lO^W/s 



(1) 



where L41 = luminosity of Cen A/10'*^erg s^^ and the 
subscript "0" refers to quantities at the source. 

Ignoring energy losses for the moment, the density of 
protons at the present time t of energy E at a, distance 
r from Cen A (assumed to be continuously emitting at a 
constant spectral rate dN^'^^ /dEdt from time ton until 
the present) is 
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where D{E) stands for the diffusion coefficient, x 
ADTon/r'^ = Ton/TD, Ton = t - ton, and 
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In each "scatter", the diffusion coefficient describes an 
independent angular deviation of particle trajectories 
whose magnitude depends on the Larmor radius i?L = 
100i;2o/S^G kpc, where E2a = E/lG'^'^eV, B^g = 
B/{1 /iG). Extragalactic magnetic field strengths and co- 
herence lengths are not well established, but it may be 
plausible to assume a Kolmogorov form for the turbulent 
magnetic field power spectrum with coherent directions 
on scales of 0.5 - 1 Mpc. One can then naively esti- 
mate that protons with energies E < lO^^^Mpc-B^G eV 
remain trapped inside magnetic subdomains of size £, at- 
taining efficient diffusion when the wave number of the 
associated Alfven wave is equal to the gyroradius of the 
particle |2|] . With a Kolmogorov spectrum this gives for 
a diffusion coefficient 
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^20 ^Mpc 
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Mpc^/Myr. 
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Here, iupc — ^/(l Mpc). For Ton — > 00, the density 
approaches its time- independent equilibrium value ricq, 
while for Ton = td — /AD, n/ricq — 0.16. 

A word about the validity of the diffusive approxi- 
mation: one may easily check that for E = 10^^ eV, 
B = O.BfiG, £ = 0.5 Mpc, the diffusive distance traveled 
ctd = 50 Mpc ^ d = 3.4 Mpc. For higher energies, the 
validity of the diffusive approach must be checked on a 
case-by case basis |Q. For these purposes, in the case of 
a continuously emitting source, the definition of a diffu- 
sion time is somewhat arbitrary. We will use td , a choice 
consistent with simulations |3l| ]. 

To further constrain the parameters of the model, 
we evaluate the energy- weighted approximately isotropic 
proton flux at 1.5 x 10^^ eV, which lies in the center of 
the flat "low energy" region |14| of the spectrum: 



{ATT)^dD{E) dEdt 
7.6 X 10^'' eL4i / eV^ m"^ g-i g^.-i 
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In the second line of the equation, we have used the val- 
ues of B and i as given in the previous paragraph. We 
fix e L41 / — 0.40, after comparing Eq.(||) to the observed 
CR-flux: E^Johs{E) = lO^^-^ eV^ m'^ s"! sr"! §. With 
eL4i ~ 1, this determines / ~ 0.40, and consequently the 
required age of the source Ton to be about 400 Myr . 
To maintain flux at the "ankle" for the same Ton, we re- 
quire an approximate doubling of Lcn at 5 x 10^^ eV. 
Because of the larger diffusive time delay at this energy, 
this translates into an increased luminosity in the early 
phase of Cen A. The associated synchrotron photons are 
emitted at energies < 30 MeV |33|. The increase in ra- 
diation luminosity in this region is not inconsistent with 
the flattening of the spectrum observed at lower energies 

In current models of describing cosmic ray acceleration, 
the principal mechanisms for energy loss are synchrotron 
radiation and photopion processes |2^,^^. If the radi- 
ation energy density of the source is sufficiently high, 
photopion production leads to copious neutron flux (that 
can readily escape the system) and associated degrada- 
tion of the proton spectrum. This occurs only near the 
maximum proton energy [ p^ . It is reasonable to assume 
that the ambient photon density of Cen A is sufficiently 
high so that near the end of the spectrum the ef- 
ficiency of neutron production e„ becomes comparable 
to the proton channel Cp. We take for granted that the 
proton spectrum cuts off at 4 x 10^" eV. Consequently, 
because of the leading particle effect Q, we expect a 
cutoff in the neutron spectrum at approximately 2 x 10^" 
eV. We adopt an energy of 1 x 10^" eV as a lower cutoff 
on the neutron spectrum, and simplify the discussion by 
assuming that in the narrow interval E20 G [1,2] e„ w ep. 
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The neutron spectrum observed at Earth is further nar- 
rowed because of decay en route. The decay length is 
A(£') = 0.9 i?20 Mpc Because of the exponential 

depletion, about 2% of the neutrons survive the trip at 
10^0 eV, and about 15% at 2 x lO^" eV. We note at this 
point that the increasing survival of neutrons at energies 
above 1.5 x 10^" eV has as a consequence of the Cen A 
model that the observed diffuse flux E'^Jobs{E) should 
begin to decrease at these energies (unless other factors 
contribute to an increase). 

We may now estimate a signal-to-noise ratio for detec- 
tion of neutron CRs in the southern hemisphere, say at 
Auger |^6|. If we assume circular pixel sizes with 2° di- 
ameters, the neutron events from Cen A will be collected 
in a pixel representing a solid angle ASl(CenA) ~ 10"'^ 
sr. For Auger {S = 3000 km^ detector with aperture 
7000 km^ sr above 10^^ eV), the event rate of (diffuse) 
protons coming from the direction of Cen A (say in a 2° 
angular cone) is found to be 



dN, 
~~dt 
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^ events/yr. 



p2 
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where we have assumed E^Jp{E) to be (approximately) 
constant up to at least E w 3x 10^° eV, in agreement with 
the observed isotropic flux in this region, E^ Jobs(-E) = 
2o24.5±o.2 eVV-^s-^sr-i 11. The neutron rate 
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events/yr, (7) 



is potentially measurable. For E20 € [1,2] we expect 



dN„ 
dt 



4 enL^i events/yr 



(8) 



arriving from the Cen A direction of the sky. With 
tnLii ~ 1/2, this gives about 2 direct events per year, 
against the negligible background of Eq.(^) |37|. Thus, 
in a few years running (of Auger) the FP hypothesis of 
Cen A as the primary source of all trans-GZK CRs can 
be directly tested. 

We now address the question of anisotropy. This can 
be found by computing the incoming current flux density 
Z?Vn as viewed by an observer on Earth, and one finds 
for a continuously-emitting source a distribution ^ (1 + 
Q!C0s(6')) about the direction of the source at angle 9 to 
the zenith, where 
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with X = Tou/td, and / was defined in Eq.(3) [Q. For 
our choices of B and £, Ton = 400 Myr, we find for E = 
eV [E = W° eV) that a = 0.04 (a = 0.07). 
It should also be remarked that the neutrons that are 
able to decay will beget secondary proton diffusion fronts 
with asymmetry parameters given by 



2D(E) 
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where 



I"{x)^ 
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du 

1,3/2 



((1 + k)u + i) 



-(1 + k)^u 



(12) 



and K = X{E)/r, X{E) being the neutron decay length 
given after Eq.(3). In spite of the complicated nature of 
Eq. ( I12I) , the results for a are very similar to the ones for 
the primary diffusion front given above. 

All in all, the Southern Auger Observatory will be in a 
gifted position to explore Cen A, providing in few years 
of operation sufficient statistics to probe extragalactic 
magnetic fields below the present observational upper 
limit 0{^G). The potential detection of the neutrons at 
Auger can subsequently be validated by the larger aper- 
ture EUSO and OWL orbiting detectors |9|. Addition- 
ally, if FP's hypothesis is confirmed, it would constitute 
a robust evidence that all FR radiogalaxies produce ex- 
tremely high energy CRs. Furthermore, our next-door 
radiogalaxy could provide a profitable arena for particle 
physics. 

In closing, we wish to comment briefly on some pub- 
lished CR observations relevant to this work. A small 
excess of fiux at 10^^ eV (detected at the Buckland Park 
field station [^) that reached the Earth preferentially 
from the direction of Cen A could militate against FP's 
hypothesis. At this energy the photon flux will be com- 
pletely damped through interactions with the cosmic mi- 
crowave background |^. Therefore, if CRs propagate 
diffusively one expects no deviation from isotropy on the 
(extragalactic) CR spectrum (except for a neutrino flux 
peaked along the line of sight). However, as far as we 
are aware, such anisotropy was not confirmed by the 
Sydney University Giant Air Shower Recorder (SUGAR) 
p2[ . Furthermore, the random arrival directions of the 
southern highest energy CRs seem to back up the above- 
outlined model. 
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